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FOREWORD 


As  part  of  the  on-going  research  program  in  "Decision 
Control  Models  in  Operations  Research",  Mr.  Carl  R.  Schultz 
investigates  the  bahavior  of  multi-item  inventory  control 
systems  at  a warehouse  level,  in  which  demand  is  comprised 
of  the  aggregated  replenishment  orders  from  lower-echelon 
inventory  control  facilities.  In  this  environment,  the 
warehouse  demand  probability  distributions  are  sporadic  and 
exhibit  correlation  from  one  time  period  to  the  next.  Under 
an  assumption  of  limited  demand  information,  Mr.  Schultz 
adapts  the  Power  Approximation  (Technical  Report  #7),  which 
was  originally  designed  for  independent  demand  distributions, 
for  warehouse  replenishment  decisions.  Several  sections  of 
this  report  parallel  similar  findings  in  earlier  reports. 
Other  related  reports  dealing  with  the  research  program  are 
listed  as  follows. 

Harvey  M.  Wagner 

Principal  Inventigator 


Richard  Ehrhardt 
Co-principal  Inventigator 
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ABSTRACT 


Inventory  managers  often  encounter  erratic  demand  histories  which 
are  difficult  to  model.  For  example,  periods  of  no  demand  are  frequently 
observed,  and  when  demand  is  positive,  it  tends  to  be  quite  large.  Fur- 
thermore, periods  of  high  demand  are  followed  by  several  periods  of  low 
demand.  One  possible  explanation  for  this  sporadic  and  correlated 
demand  behavior  is  that  demand  originates  from  separate  facilities  which 
employ  (s,S)  replenishment  policies.  That  is,  each  period's  demand 
is  the  sum  of  the  replenishment  order  quantities  received  from  other 
facilities.  In  this  investigation  we  consider  the  management  of  an 
inventory  system  exposed  to  this  type  of  demand  environment. 

The  complex  nature  of  the  demand  history  that  arises  in  such  a 
wholesale  warehouse  environment  makes  the  computation  of  optimal  replen- 
ishment policies  prohibitive.  We  confine  our  attention,  therefore,  to 
stationary  policies  of  the  (s,S)  form,  and  estimate  the  best  values 
of  the  policy  parameters  by  means  of  computer  simulation. 

In  an  applied  setting  it  is  impractical  to  use  simulation  to 
derive  an  inventory  policy  for  each  item  of  interest.  We  therefore  seek 
an  approximately  optimal  policy  rule  that  is  easy  to  compute.  Two  policy 
rules  are  examined  in  the  wholesale  warehouse  environment.  In  order  to 
discover  the  cost  implications  of  ignoring  the  sporadic  and  correlated 
nature  of  the  warehouse  demand  process,  we  first  examine  a stationary 
(s,S)  policy  rule  that  does  not  recognize  these  demand  properties.  For 
our  experimental  design,  the  policy  rule's  total  cost  performance  is 
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unsatisfactory;  typically,  8%  above  that  of  stationary  (s,S)  policies 
found  by  simulation.  The  other  policy  rule  we  examine  is  an  adaptation 
of  the  Power  Approximation  of  Ehrhardt  (1976)  to  a correlated  demand 
environment.  The  only  demand  information  required  by  the  policy  rule 
is  the  mean,  variance,  and  variance  over  one  leadtime  of  demand.  For 
the  same  experimental  design,  total  cost  performance  is  typically  only 
above  that  of  policies  found  by  simulation. 
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1.  (s,S)  INVENTORY  POLICIES  FOR  A WHOLESALE  WAREHOUSE  SYSTEM 


In  many  inventory  systems  one  often  observes  an  erratic  demand 
history  whose  underlying  generating  process  is  difficult  to  explain. 

For  example,  consider  a situation  in  which  periods  of  no  demand  are 
frequently  observed,  and  when  demand  is  positive,  it  tends  to  be  quite 
large.  One  possible  explanation  for  this  behavior  is  that  demand  origi- 
nates from  other  facilities  which  employ  (s,S)  replenishment  policy 
rules.  If  this  is  true,  demand  may  also  be  correlated  from  one  period 
to  another.  We  consider  the  use  of  (s,S)  replenishment  policies  for  an 
erratic,  correlated  demand  history  that  arises  in  a wholesale  warehouse 
environment.  We  analyze  a two-echelon  inventory  system  consisting  of 
a number  of  lower-echelon  facilities  (stores)  satisfying  exogenous  cus- 
tomer demand  and  themselves  acting  as  customers  to  a single  upper-echelon 
facility  (warehouse). 

Most,  if  not  all,  of  the  research  associated  with  the  stochastic, 
two-echelon  inventory  control  problem  [Clark  (1972)],  has  been  in  the 
direction  of  establishing  rules  or  policies  which,  when  applied  at  the 
individual  facilities,  satisfy  a prescribed  system-wide  performance 
objective,  such  as  minimal  expected  cost.  The  complexity  of  centralized 
system  control  makes  the  problem  of  determining  optimal  order  and  supply 
policies  computationally  untractable  unless  very  restrictive  assumptions 
are  imposed  on  both  the  model  and  the  policy  forms,  assumptions  which 
tend  to  be  over-restrictive  relative  to  practical  applications.  In  addi- 
tion, centralized  control  over  the  entire  system  is  often  not  possible 
in  applied  situations.  In  this  investigation  we  depart  from  the  tradi- 
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tional  approach  of  system-wide  inventory  control  and  instead  view  the 
management  of  the  warehouse  simply  as  a single-facility  control  problem. 
The  demand  process  at  the  warehouse  is  given  by  the  aggregated  replen- 
ishment processes  of  the  stores  which,  we  shall  assume,  employ  (s,S) 
policies.  As  noted  earlier,  this  demand  process  can  be  highly  erratic 
and  significantly  correlated,  unlike  the  demand  processes  commonly 
studies  in  the  single-facility  inventory  control  problem.  In  fact, 
little  research  has  dealt  with  the  case  of  dependent  demands,  and 
existing  results  [Karlin  and  Fabens  (1959)]  are  not  applicable  to  this 
demand  structure. 

1 . 1 A Wholesale  Warehouse  Inventory  Model 

We  consider  the  management  of  a wholesale  warehouse  inventorv 
system.  We  model  the  wholesale  warehouse  as  the  upper-echelon  facility 
of  the  two-echelon  inventory  system  graphically  represented  in  Figure  1.1 

Basic  demand,  which  we  assume  is  independent  and  identically  dis- 
tributed, enters  the  two-echelon  inventory  system  at  a lower-echelon 
comprised  of  N independently-operated  stores  in  parallel.  The  basic 
demand  is  filtered  through  the  replenishment  policies  at  each  store  and, 
disallowing  transshipments  of  items  between  the  stores,  is  passed  to  the 
warehouse  in  the  form  of  aggregated  replenishment  orders.  Specifically, 
if  X.j  represents  the  order  quantity  received  from  store  k in  period 
i , then 

N k 

(1.1)  Zi  = I Xi 

1 k=l  1 

is  the  entire  demand  realized  at  the  warehouse  in  period  i . 


We  postulate  a single-item,  periodic  review  inventory  model  for 
each  facility  in  the  two-echelon  system.  We  assume  that  demand  at  each 
facility  is  met  as  long  as  there  is  stock  on  hand,  and  when  a stockout 
occurs,  unfilled  demand  is  backlogged  until  sufficient  replenishments 
arrive.  Items  kept  in  inventory  are  assumed  to  be  conserved,  there 
being  no  losses  by  deterioration,  obsolescence,  or  pilferage.  Inventory 
on  hand  at  the  end  of  a current  period  is  the  inventory  from  the  pre- 
vious period  plus  any  replenishment  that  arrives,  less  demand  in  the 
current  period.  If  inventory  on  hand  is  negative,  its  absolute  value 
is  the  amount  of  backlogged  demand.  The  time  sequence  of  events  within 
any  period  is  taken  to  be  order,  delivery,  and  demand. 

The  cost  structure  at  the  warehouse  is  of  a simple  form.  At  the 
end  of  each  period  a unit  holding  cost  h or  a unit  penalty  cost  p 
is  incurred  for  each  unit  on  hand  or  on  back  order,  respectively.  The 
cost  of  a replenishment  quantity  is  assumed  to  be  linear  with  a fixed 
setup  cost  K , and  warehouse  replenishments  are  assumed  to  be  delivered 
a fixed  lead  time  L after  being  ordered.  The  criterion  for  optimal 
inventory  control  is  minimization  of  the  undiscounted  expected  cost  per 
period  over  an  infinite  horizon. 

We  assume  that  control  over  replenishment  at  each  facility  in  the 
two-echelon  system  is  exercised  by  an  (s,S)  policy:  whenever  inventory 
x on  hand  and  on  order  is  less  than  or  equal  to  the  value  s , an  order 
is  placed  for  a replenishment  of  size  S - x . 

Under  the  model  assumptions  described  above,  an  (s,S)  policy  is 
known  to  be  optimal  if  demands  are  independent  and  identically  distributed 


id 


, 
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[Iglehart  (1963a,  b)]  , and  methods  for  computing  optimal  policies 
have  been  developed  [Veinott  and  Wagner  (1965)]  . When  demands  in 
successive  periods  are  dependent,  as  is  the  case  at  the  warehouse,  an 
optimal  policy  will  not  be  of  the  (s,S)  form.  Nevertheless,  in  this 
study  we  confine  our  attention  to  warehouse  replenishment  policies  of 
the  (s,S)  form,  since  their  simple  form  has  led  to  their  frequent 
use  in  applied  situations. 

1 .2  An  Algorithm  to  Set  (s,S)  With  Partial  Knowledge 

We  base  our  algorithm  to  determine  values  for  the  warehouse  policy 
parameters  s and  S on  the  Power  Approximation  of  Ehrhardt  (1976)  . 

The  Power  Approximation  is  an  algorithm  for  computing  approximately 
optimal  values  for  (s,S)  using  only  the  mean  y and  variance  o2  of 
demand  and  is  based  on  asymptotic  renewal  theory.  Furthermore,  the  Power 
Approximation  is  derived  for  independent,  identically  distributed  demands. 
One  of  the  goals  of  this  study  is  to  investigate  an  adaptation  of  the 
Power  Approximation  to  the  sporadic  and  correlated  demand  structure 
present  in  our  wholesale  warehouse  inventory  system. 

The  Power  Approximation  is  executed  as  follows.  Let 

(1.2)  Dp  = (1.463)y-364(K/h)'498(oL+1)-1382 
and 

(1.3)  sp  = (L  + l)y  + (aL+1)’832(a2/u)J87(.220/z  + 1.142  - 2.866z)  , 


where 
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(1.4)  z = {Dp/C ( 1 + p/h)aL+1]}'5  , 
and 

(1.5)  o^+l  * (L  + 1)fl2  • 

If  Dp/p  is  greater  than  1.5  , let  s = sp  and  S = sp  + Dp  . Other- 
wise compute 

(1.6)  SQ  = (L  + l)u  + voL+1  , 
where  v is  the  solution  to 

(1.7)  f ( 2tt )-1s  exp(-x2/2)  dx  = p/(p  + h)  . 

-oo 

The  policy  parameters  are  given  by 

s = minimum(Sp,  Sq) 

(1.8) 

S = minimum(Sp,  SQ) 

If  demands  are  integer-valued,  then  sp  , Dp  , and  Sq  are  rounded 
to  the  nearest  integer. 

Since  the  Power  Approximation  is  designed  for  independent,  identi- 
cally  distributed  demands,  the  term  a£+-|  given  by  (1.5)  represents 
the  variance  of  demand  over  L+l  periods.  Our  adaptation  of  the  Power 
Approximation  adjusts  this  term  to  account  for  correlation  in  the  demand 


process. 
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If  the  demand  process  is  covariance-stationary,  then  the  variance 

, 2 

aL+i  of  demand  over  L+l  periods  can  be  written  as 

(1.9)  o*  = [(L  + 1)  + 2 l (L+l  - j)p(j)]a2  , 

L+l  j=l 

where  p(j)  is  the  correlation  between  demands  separated  by  j periods. 
Thus  our  adaptation  of  the  power  approximation  is  the  replacement  of 
(1.5)  with  (1.9)  . 

In  Section  2 we  present  a detailed  analysis  of  the  demand  process 
at  the  warehouse.  Expressions  for  the  mean,  variance,  and  autocorrela- 
tion function  of  the  warehouse  demand  process,  in  terms  of  the  store 
parameters,  are  derived. 

In  Section  3 we  discuss  several  cost-effective  policy  rules  for 
the  wholesale  warehouse.  The  complexity  of  the  warehouse  demand  process 
makes  the  computation  of  optimal  policies  prohibitive.  We  confine  our 
attention  to  policies  of  the  (s,S)  form  and  use  simulation  to  estimate 
the  minimum-cost  (s,S)  policy  for  use  in  the  warehouse  environment. 

The  simulation  program  which  finds  these  "best"  (s,S)  policies  is  also 
described  in  Section  3 . Finally,  an  experimental  design  is  presented 
for  the  evaluation,  in  later  sections,  of  approximately-optimal  policy 
rules.  The  performance  of  these  policy  rules  in  the  wholesale  warehouse 
environment  is  evaluated  by  comparing  their  key  operating  characteristics 
with  the  corresponding  operating  characteristics  of  best  (s,S)  policies 
under  the  same  conditions. 

In  Section  4 an  evaluation  of  policy  rules  that  ignore  the  sporadic 
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and  correlated  properties  of  the  warehouse  demand  process  is  presented. 
In  particular,  these  are  policy  rules  that  would  be  optimal  if  warehouse 
demands  were  independently  and  identically  distributed  with  a negative 
binomial  distribution.  We  find  the  performance  of  these  policies  to 
be  unsatisfactory.  Typically,  for  a large  multi-item  system,  expected 
total  costs  are  approximately  8%  higher  than  those  associated  with 
the  simulation-derived  best  (s,S)  policies. 

In  Section  5 an  evaluation  of  the  Adjusted  Power  Approximation, 
discussed  above,  is  presented.  We  find  that  by  adjusting  for  correla- 
tion, significant  improvement  can  result.  Specifically,  expected  total 
costs  are  about  3%  above  best  values  for  the  same  multi -item  system. 

Section  6 concludes  this  investigation  with  a summary  of  our  find- 
ings, conclusions,  and  suggestions  for  future  research. 
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2.  ANALYSIS  OF  THE  WAREHOUSE  DEMAND  PROCESS 

In  this  section  we  present  a detailed  analysis  of  the  demand 
environment  in  the  wholesale  warehouse  inventory  model  described  in 
Section  1. 

In  Subsection  2.1  we  study  the  replenishment-quantity  process 
of  a single  store  employing  a stationary  (s,S)  policy  and  experienc- 
ing independent,  identically  distributed,  integer-valued,  non-negative 
demands.  In  particular,  we  develop  the  stationary  distribution  of  the 
process,  and  derive  simple  expressions  for  its  mean  and  variance.  We 
also  present  a recursive  formula  for  the  autocorrelation  function  of 
the  process. 

Finally,  in  Subsection  2.2  we  make  use  of  our  model  assumptions 
to  derive  simple  expressions  for  the  mean,  variance  and  autocorrelation 
function  of  the  warehouse  demand  process  in  terms  of  their  counterparts 
for  the  single  store  replenishment-quantity  process. 

2.1  Properties  of  the  Replenishment-Quantity  Process  at  a Single  Store 

We  assume  the  demand  process  at  a store,  denoted  by  q^ , q£,  •••  , 
is  a sequence  of  non-negative,  integer-valued,  independent  and  identically 
distributed  random  variables  having  cumulative  distribution  $(•)  and 
probability  mass  function  $(•)  . Let  4>n(-)  and  <j>n(*)  be,  for  n > 1 , 
the  n-fold  convolutions  of  $>(•)  and  $>(•)  . Let  $>^(*)  represent 


the  distribution  whose  full  mass  is  located  at  zero.  That  is, 
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(2.1) 


which  implies  that 


(2.2) 


4>°(y) 


<t>°(y) 


We  define  the  renewal  functions 


(2.3) 


and 


M(y)  = l <1>  (y)  , 

k=l 


y < 0 
y ^ 0 , 


y * 0 
y = 0 . 


(2.4) 


m(y)  = l <t>k(y) 

k=l 


1 


Let  us  and  denote,  respectively,  the  mean  and  variance  of 
the  store's  demand  distribution. 

We  assume  that  the  store  employs  a stationary  (s,S)  replenishment 
policy.  Let  D = S - s and  note  that  the  possible  replenishment  quantity 
values  are  0,  D+l , D+2,  ...  . 

Using  the  store  assumptions  given  above,  we  analyze  the  properties 
of  the  store's  replenishment-quantity  process.  Several  important  results 
from  this  analysis  are  given  by  the  following  lemma. 
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Lemma  2.1: 

If  X.  represents  the  replenishment  order  quantity  in  period  i , 

then 

i)  the  stationary  distribution  of  X.  is  given  by 
Pr[X.  = 0]  = M(D)/[1  + M(D)]  , 
and  for  k = 0,1 , 2,  . . . , 

D 

Pr[X.  = D+k+1  ] = [<j>(D+k+l ) + \ <j>(  j+k+1  )m(D-j )]/[!  + M(D)]  , 

j=0 

ii)  the  expected  value  of  X,.  , denoted  by  , satisfies 


ur  = ws 


iii)  the  variance  of  X.  denoted  by  , satisfies 
? ? ^ 

of  = + [2p_  l km(k)]/[l  + M(D)]  , 

r s s k=C 


and 


M 


Proof: 

Let  represent  the  replenishment  order  quantity  in  period 
i . We  proceed  to  derive  an  expression  for  the  distribution  of  Xi 
conditioned  on  the  number  of  periods  n since  the  last  order  was 
placed. 

First,  consider  the  case  when  an  order  was  placed  in  the  preced- 
ing period,  that  is,  n = 1 . Then  for  0 < t < D 


(2.5) 


Pr[Xi  s 1 1 n = 1]  = Pr[Xi  = 0 1 n = 1] 
= Pr[q1  < D] 

= *(D)  , 


and  for  t > D , 


(2.6) 


Pr[Xi  < t | n = 1]  = Pr[q]  < t] 

= *(t)  . 

Now  suppose  n > 1 . For  0 < t s D we  have 

(2.7)  Pr[X.  < t | n = j]  = Pr[X.  = 0|n  = j] 

= Pr[q]+q2+. . .+q^  s Dlq^q^. . .+qj._]  <;  D] 
= *j(D)/*j_1(D)  . 

For  t > D , we  have 

(2.8)  Pr[X.  s t | n = j]  = Pr[X^  s D|n  = j]  + Pr[D  < Xi  < t|n=j] 

= (4>J(D)  + Pr[D  < q^+q2+. . .+qj  s t 

and  q1+q2+...+qj_1  s D]}/4>J-1  (D) 
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Observe  that 


(2.9)  Pr[D  < q1+q2+...+qj  < t and  q1+q2+. . .+qj_1  s D] 

= Pr[q-|+q2+.  • - +q j _ -j  * min(D,t  - q^)]  - Pr[q1+q2+. . •+qj_1  * D - q^] 
Conditioning  on  the  value  of  q.  , the  right-hand  side  of  (2.9)  becomes 

J 

(2.10)  *J’"1(D)*(t-D-l)  + [ (t-k)4>(k)  - 4>j(D)  . 

k=t-D 


Therefore, 


(2.11)  Pr[Xi  <;  t | n = j]  = $(t-D-l ) + 


l ( t-k)4>(  k) 

k=t-D 


/*j"1(D) 


To  summarize,  for  j = 1,2,3,..., 


(2.12) 


Pr[Xi  < t|n  = j] 


0 

^(D)/<tJ_1 
4>  ( t-D- 1 ) + 


t < 0 
0 < t < D 
( t-k  )d>(  k )! (D) 


t > D 


Using  (2.12)  we  can  readily  obtain  the  following  conditional  proba- 
bilities: 

Pr[X.  = 0 1 n = j]  = 4>j(D)/4>j_1(D)  , 


(2.13) 


j = 1 ,2,3,. . . 
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1 


1 t 


(2.14)  Pr[X.  = D+m+1 | n = j] 


= l 4>(k+m+l)4>j‘^D-k)j/4>j 


•1 


(D)  , 

m = 0,1,2,... 


Using  a Markov  chain  analysis,  we  can  determine  the  stationary  or 
limiting  distribution  of  the  number  of  periods  n between  orders. 
Specifical ly. 


(2.15) 


Pr[n  = j]  = *j“1(D)/[l  + M(D)]  , j = 1,2,3,... 


The  limiting  probabilities  given  in  (2.15)  enable  us  to  explicitly 
determine  the  stationary  distribution  of  the  single  period  replenishment 
quantity  X.  . For  m = 0,1,2,...  , we  have 


(2.16)  Pr[X . = D+m+1]  = l Pr[X.  = D+m+1 |n  = j]  Pr[n  = j]  . 
1 j=l  1 


Using  (2.14)  we  can  express  the  right-hand  side  of  (2.16)  as 


(2.17) 


D . , 

l l 4>(  k+m+1  )<jr  (D-k) 
_j  = l k=0 


<j)(  D+m+1 ) + l 4>(  k+m+1 ) l <j>J(D-k) 
k=0  j=l 


/[l  + M(D)  ] 

/[l  + M(D)]  . 


Thus 


(2.18)  Pr[Xi  = D+k+1]  = 


D 


d>  ( D+k+1 ) + l <(>(  j+k+1  )m(D-j ) 
j=0 


/[l  + M(D)]  , 


k = 0,1,2,... 


In  a similar  fashion  we  obtain 


(2.19) 


Pr[X.  = 0]  = M(D)/[1  + M(D)] 


15 


i 

P 

i 

I 

I 
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Thus  part  i has  been  established. 

Since  all  demand  eventually  finds  its  way  into  replenishment  orders, 
we  would  expect  the  mean  of  the  stationary  replenishment  distribu- 
tion, given  by  (2.18)  and  (2.19),  to  be  the  same  as  the  mean  us  of 
the  store's  demand  distribution.  This  intuitive  result  can  be  shown 
directly  using  the  stationary  replenishment  distribution.  The  proof, 
however,  is  not  presented  here. 

As  part  ii  of  the  lemma  indicates,  the  mean  ur  of  the  stationary 
replenishment  distribution  does  not  depend  on  the  store's  policy  param- 
eters s and  S . As  one  might  well  anticipate,  however,  the  variance 
2 

a of  the  stationary  replenishment  distribution  depends  upon  D . 

Let  X.j  represent  the  single  period  replenishment  quantity.  Then 
from  part  ii  of  the  lemma 

(2.20)  Var  X.  = EX?  - (EX.)2 


We,  therefore,  need  to  evaluate 


(2.21) 


? 00  o 

EXT  = l (D+n+1 } Pr[X.  = D+n+1 ] 
1 n=0  1 


Using  (2.18)  we  have 


(2.22)  EXf  = 


9 u 00  ? 

I (D+n+1  )^<|>(D+n+l ) + l m(D-y)  £ (D+n+1 ) 4>(y+n+l ) 


n=0 


y=0 


n=0 


[1  + M ( D ) ] 
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Observe  that 


(2.23)  l (D+n+l)2<f>(D+n+l)  = a2  + - l k24>(k)  , 

n=0  k=0 


where  o2  is  the  variance  of  the  store's  demand  distribution, 
s 


(2.24)  l (D+n+1  )24>(y+n+l ) = l [(D-yr  + 2(D-y) (y+n+1 ) + 
n=0  n=0 


= (D-y)2  + °s  + us  + ' 


With  the  aid  of  (2.23)  , (2.24)  , and  the  fact  that 


(2.25) 

we  can  rewrite  (2.22)  as 


M(z)  = l m(j)  , 
3=0 


(2.26)  EX^  ' “s  * "s  + [2"s j0  km(k)]/['  + H<D)] 

(o  o y 9 o ? ) 

+ ] l m(D-y)[(D-y)2  - l (y-D-k)^(k)]  - l k‘<t>(k) 
(y=0  k=0  k=0  ) 


Rearranging  terms,  we  have 

l m(D-y)[(D-y)2  - f (y-D- k ) 2<f. ( k ) ] - l k2<J>(k) 


(2.27) 


y=o 


k=0 


k=0 


= l (D-y)2[m(D-y)  - l ( k )m( D-y-k ) 
y=0  k=0 


Also 

(y+n+1 )2] 

4>  (y+n+1 ) 

f (y-D-k)2 
k=0 

4>(  k ) 


/[I  + M(D)1  . 


- 4>(D-y )] 


which  is  zero,  since  m(*)  satisfies  the  renewal  equations 


(2.28) 


Therefore, 


(2.29) 


m(z)  = <p(z)  + l c; ' v)m(z-y) 
y=0 


for  z > 0 


+ [2u  l km(k)]/[l  + M(D)]  , 

s s k=0 


and  part  iii  of  the  lemma  has  been  established.  It  can  also  be  shown, 

2 

as  one  would  expect,  that  or  is  increasing  in  D . 

We  previously  observed  that  due  to  the  nature  of  (s,S)  policies 
the  replenishment  process  {X^;  i > 1}  will  exhibit  autocorrelated 
behavior.  Since  the  replenishment  process  is  covariance-stationary, 
the  covariance  between  X.  and  X.+.  will  be  independent  of  i . The 
autocovariance  of  the  replenishment  process  at  lag  j , denoted  y (j)  , 
is  defined  as 

(2.30)  yr(j)  = Cov(X.,  X.+j)  = E[(X.  - ur)(Xi+j  - Mr)] 

= Exixi*j  - 4 • 

We  have 

(2.31)  EX.X.+.  = E[X.Xi+.|X.  = 0]  Pr[X.  = 0]  + E[X.X.+j|X.  > 0] 

• Pr[Xi  > 0]  . 

Notice  that 


(2.32) 


ECVhjIX,  • 0]  ■ 0 . 
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Furthermore,  each  positive  X,  marks  a regeneration  point  in  the 

J 

process.  Thus  if  X.  > 0 , X^  and  X^+j  will  be  independent  for 
j = 1,2,3,...  . Consequently, 

(2.33)  E[X.X.+j|X.  > 0]  = E[X.|X.  > 0]  E[X.+j|Xi  > 0] 

Using  (2.32)  and  (2.33)  , (2.31)  can  be  rewritten  as 

(2.34)  EX.X.+.  = E[X.|X.  > 0]  Pr[X.  > 0]  E[Xi+j|X.  > 0] 

=UrE[Xjtj|X(  >0]  . 

Therefore, 

(2.35)  Yr(j)  = ur{E[Xi+j|Xi  > 0]  - ur}  . 

To  evaluate  E[X.+.|X.  > 0]  , we  need  an  expression  for  the  condi- 
tional  expectation  of  X^  given  that  n periods  have  elapsed  since  an 
order  was  last  placed.  Using  the  conditional  probabilities  given  in 
(2.13)  and  (2.14)  , we  have 


(2.36) 


E[X . | n = j]  = l (D+m+1 ) Pr[X.  = D+m+l|n  = j] 
1 m=0  1 


D i 1 

l f3  ' (D-y)  I 
y=0  m=0 


(D+m+1  )4>(y+m+l ) 


/<Kj_1(D) 


Observe  that 


(2.37)  J (D+m+1  )4>(y+m+l ) 
m=0 


y 

v.  + D - y + I (y-D-k)<(>(k) 
s k=0 
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Rearranging  terms  yields 


(2.38)  l (J)'’-1  (D-y)  l (y-D-k)<J>(k)  = - l (D-y)  l <J>(k)<^_1  (D-y-k) 
y=0  k=0  y=0  k=0 


Note  that  the  right-hand  side  of  (2.38)  is  equivalent  to 


(2.39) 


■ I M>J(k)  . 
k=0 


Thus  we  have 


(2.40)  E[X.|n  = j]  = ur  + 


\l 

Lk=0 


-1 


(k)  - <J>j(k)) 


/*j_1(D) 


For  j 
(2.41) 


= 1 


E[Xi+1|Xi  > 0]  = E[X.+1|n  = 1]  = wr  - l k<fr(k) 


For  j = 2 


(2.42)  E[X.+2|Xi  > 0]  = E[X.+2|X.  > 0 and  X.+]  = 0]  Pr[X.+1  = 0|Xi  > 0] 


+ E[X.+2|X.  > 0 and  X.+]  > 0]  Pr[Xi+]  > 0|X.  > 0] 


or  equivalently 


(2.43)  E[Xi+2|X.  > 0]  = E[X.+2|n  = 2]  Pr[X.+]  = 0|n  = 1] 


+ E[X.+2|n  = 1]  Pr[Xi+1  > 0|n  = 1] 


Employing  expressions  (2.13)  , (2.40)  , and  (2.41)  we  obtain 


j 


y 
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(2.44)  E[X..JX.  > 0]  = m + 4(D)  £ M>(k)  - £ k^(k) 

1 * 1 r k=0  k=0 


= Mr  - ][  k<t>  (k)  + 4(D){yr  - E[x -+1 1 X - > 0]} 


Continuing  in  this  fashion,  we  have  the  following  recursive  formrla 
for  E[X.+j|X.  > 0]  : for  j = 1,2,3,...  , 


D 


ur  - I M>(k) 
r k=0 


j = 1 


(2.45)  E[X.+j|X.  > 0]  = J 


wr  - l k4>J'(k)  + V ^(D) 


k=0  i= 1 

K - ECXi+e|Xi  > 0]} 


j > 1 


By  employing  (2.35)  we  have 


(2.46) 


Yr(j) 


= < 


-ur  l k4>(k) 


j = 1 


k=0 


-ur  l k4J(k) 
r k=0 


j-1  H t 

l 4J  *(D)v  U) 
1=  1 r 


j > 1 


If  we  let  p (j)  denote  the  autocorrelation  at  lag  j of  the 


replenishment  process  (X.;  i > 1}  , then  substitution  of 


2 

Pr(j)  = Yr(j)/°r 


(2.47) 
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into  (2.46)  yields  part  iv  of  the  lemma,  and  the  proof  is  complete. 

Notice  that  the  first-order  autocorrelation  Pr(l)  is  negative 
and  becomes  more  negative  as  D increases.  Both  of  these  observations 
are  not  unexpected  considering  the  nature  of  (s,S)  policies.  A 
replenishment  order  in  one  period  is  followed  by  several  periods  with 
low  probability  of  orders  being  placed,  and  as  D becomes  larger, 
replenishment  orders  become  less  likely. 

2.2  Properties  of  the  Warehouse  Demand  Process 

In  the  previous  subsection  we  developed  expressions  for  the  mean 
2 

and  variance  or  of  the  stationary  replenishment  process  and  for 
the  correlation  between  replenishment  orders  separated  by  j periods, 
Pr(j)  • The  properties  of  our  wholesale  warehouse  inventory  model, 
described  in  Section  1 , enable  us  to  derive  simple  expressions  for 
the  mean,  variance  and  autocorrelation  function  of  the  warehouse  demand 
process  in  terms  of  their  counterparts  for  the  single-store  replenish- 
ment processes. 

o 

Let  ur>k  , ar  k , and  yr  k(-)  be,  respectively,  the  mean, 

variance  and  autocovariance  function  of  the  replenishment  process  at 

2 

store  k . Let  y , a , and  p (•)  be,  respectively,  the  mean, 
w w w 

variance,  and  autocorrelation  function  of  demand  process  at  the  warehouse. 
Several  important  properties  of  the  warehouse  demand  process  are  given 
in  the  following  lemma. 

Lemma  2.2: 

A.  If  Z..  represents  the  demand  at  the  warehouse  in  period  i. 


1 


1 * 

i}  ww  = ^ ur,k 

2 v 2 

”>  “w'j,  ’’r.k 

N N 2 

iii)  pw(j)  - ^ Yr>k(J)/k|i  °r,k 


for  j = 1,2,3,...  ; 


B.  Furthermore,  if  each  store  has  the  same  demand  distribution 
and  replenishment  policy,  then 


iv)  uw  = Nyr 

v)  o*  = Na^ 


vi)  pw(j)  = Pr(j)  . 


for  j = 1 ,2,3, . 


Proof: 


Let  Zi  represent  the  demand  at  the  warehouse  in  period  i 


(2.48) 


k 

Zi  = l x-  . 

1 k=i  1 


where  is  the  order  quantity  received  from  store  k in  period  i 


Thus, 


(2.49) 


ez'  • il x']  ■ 


I EX  , 
k=l  1 


and  part  i is  established.  Since  the  stores  operate  independently  of 


each  other  we  have 
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(2.50) 


Var  Z.  = Var 


U x<]  ■ j, 


Var 


and  part  ii  is  established. 


Let  Yw(j)  denote  the  autocovariance  at  lag  j of  the  warehouse 


demand  process,  that  is,  Yw(j)  is  defined  as 


(2.51) 


(j)  s Cov(Z . , 1M)  = EZ .Z . . . - 


’i*  i+j 


1 i+j  Mw 


From  (2.48)  and  (2.49)  we  have 


(2-52)  x»(j)  = {j, j,  < J - [ J,  v J 


N k k N N k / 

/KyK  \ • r r r/vKy-C 


l E(x;xj+1)  + l l E(XjX*  ) - l V k 
<=1  1 1 J k=l  1=  1 1 1 J k=l  r,K 


k * t 


N N 

l l v .w  / 
k=l  £=1  r’K  r,t 
k ^ l 


The  independent  operation  of  the  stores  allows  us  to  rewrite 
(2.52)  as 


<2-53>  ■ j,  tE<x?xHj)  - V.kl  ■ j,  Vk(j>  • 


Thus  we  have  part  iii  of  the  lemma. 

Note  that  if  each  store  has  the  same  demand  distribution  and  replen- 
ishment policy,  then  the  replenishment  processes  will  be  identically 
distributed  and  part  B of  the  lemma  follows  directly. 


* i 
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3.  EVALUATING  THE  PERFORMANCE  OF  WAREHOUSE  INVENTORY  POLICIES 

In  this  section  we  discuss  several  cost-effective  policy  rules 
for  the  wholesale  warehouse  inventory  model  and  describe  the  methods 
used  to  evaluate  their  performance. 

3.1  Warehouse  Policy  Rules 

Table  3.1  lists  several  cost-effective  replenishment  policies 
for  a wholesale  warehouse  environment.  Brief  comments  about  the 
nature  of  the  policy  rules  and  the  sections  in  which  they  are  studied 
are  also  given. 

Recall  that  our  criterion  for  optimal  warehouse  inventory  control 
is  minimization  of  the  undiscounted  expected  cost  per  period  over  an 
infinite  horizon.  The  complexity  of  the  warehouse  demand  process  makes 
the  computation  of  that  policy  which  meets  this  criterion  prohibitive. 
Even  the  form  of  the  optimal  policy  is  very  difficult  to  characterize. 

The  optimal  policy  will  not  be  of  a stationary  (s,S)  form.  Neverthe- 
less, we  have  confined  our  attention  to  the  (s,S)  form  because  of  its 
popular  use. 

Even  if  we  restrict  the  cost  minimization  search  to  only  those 
policies  of  a stationary  (s,S)  form,  the  computational  difficulties 
remain.  Hence,  we  investigate  approximately  optimal  (s,S)  policy  rules 
Any  such  investigation  of  approximately  optimal  (s,S)  policy  rules, 
however,  will  be  fruitless  unless  there  is  a benchmark  by  which  we  can 
evaluate  their  performance.  Ideally,  this  benchmark  is  provided  by  the 
operating  characteristic  values  of  an  optimal  stationary  (s,S)  policy. 


Table  3.1 

Replenishment  Policies  For  A Wholesale  Warehouse  Environment 


POLICY  RULE 

COMMENTS 

OPTIMAL 

Complex  to  characterize 

(not  of  a stationary  (s,S)  form). 

MINIMUM  COST 
STATIONARY 
(s,S) 

Complex  to  compute. 

VERY  BEST 
(s,S) 

Minimum-cost  stationary  (s,$) 
policy  for  a long  finite  his- 
tory of  demand. 

BEST  (s ,S) 

Fibonacci  search  for  very  best 
(s,S)  policy. 

I.I.D.  OPTIMAL 

Optimal  stationary  (s,S)  policy 
for  an  i.i.d.  demand  process  fol- 
lowing a negative  binomial  distri- 
bution having  the  same  mean  and 
variance  as  the  warehouse  demands. 

APPROXIMATELY 

OPTIMAL 

Power  Approximation  corrected  for 
correlation. 

SECTION 
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Unfortunately,  these  operating  characteristic  values  also  are  prohibi- 
tively difficult  to  compute. 

To  circumvent  all  these  difficulties  we  estimate  an  optimal  sta- 
tionary (s,S)  policy  by  simulation.  We  seek  the  stationary  (s,S) 
policy  that  minimizes  the  total  cost  per  period  for  a long  history  of 
generated  warehouse  demands;  we  refer  to  this  policy  as  the  "very  best" 
(s,S)  policy.  The  software  package  used  to  find  very  best  (s.S) 
policies  was  developed  by  Kaufman  (1976)  . The  methodology  of  his 
program  is  outlined  in  the  following  paragraphs. 

A sequence  of  warehouse  demands  is  generated.  A sequence  of 
warehouse  (s.S)  policies  is  empirically  tested  through  simulation  on 
the  generated  demand  sequence,  and  a policy  with  the  smallest  total 
cost  per  period  is  selected  as  the  very  best  (s,S)  policy. 

More  specifically,  a sequence  of  values  for  D is  examined.  For 
each  D , a corresponding  value  for  S is  selected  that  minimizes  the 
total  cost  per  period  given  the  generated  demand  sequence.  The  sequence 
of  D values  is  determined  by  a Fibonacci  search  over  a range  of  possi- 
ble values  for  D [Wagner  (1968)]  . 

This  technique  would  guarantee  finding  the  very  best  (s.S)  policy 
for  the  generated  demand  sequence  if  the  total  cost  per  period  were 
convex  in  D ; but  such  convexity  is  generally  not  present  [Wagner, 
O'Hagan,  and  Lundh  (1965)]  . Thus,  the  policy  found  may  not  be  the 
very  best  (s.S)  policy  for  the  generated  demand  sequence.  We  refer 
to  the  policy  resulting  from  the  search  technique  described  above  as 
the  "best"  (s,S)  policy. 
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Once  a best  (s,S)  policy  has  been  found,  we  operate  the  policy 
on  the  same  generated  demand  history  and  collect  key  operating  charac- 
teristic values.  Since  the  computed  operating  characteristic  values 
are  statistical  estimates,  confidence  intervals  are  also  computed. 
Kaufman's  program  also  has  an  option  that  permits  the  computation  of 
operating  characteristic  values  and  confidence  intervals  for  an  arbi- 
trary (s,S)  warehouse  policy.  By  using  this  option  we  can  evaluate 
the  performance  of  other  (s,S)  policies  by  comparing  their  operating 
characteristic  values  against  those  of  the  best  (s,S)  policy  on  an 
identical  warehouse  demand  sequence. 

Although  the  simulation  program  discussed  above  computes  nearly 
optimal  (s,S)  policies,  its  implementation  in  an  actual  wholesale  ware- 
house inventory  system  is  impractical.  It  requires  the  knowledge  of 
the  demand  distribution  and  policy  parameters  at  each  store,  informa- 
tion usually  unavailable  to  the  warehouse  manager.  In  addition,  for 
a large  multi-item  system  this  method  could  be  very  expensive.  For 
these  reasons  we  seek  policy  rules  that  require  only  limited  warehouse 
demand  information. 

The  first  rule  of  this  type  that  we  consider  is  one  which  ignores 
the  sporadic  and  correlated  nature  of  the  warehouse  demand  process. 

This  policy  rule  would  be  optimal  at  a facility  having  the  same  cost 
structure  as  the  warehouse,  but  with  an  independent  and  identically 
distributed  demand  process  following  a negative  binomial  distribution 
having  the  same  mean  and  variance  as  the  warehouse  demands.  We  refer 
to  this  policy  rule  as  the  I.I.D.  optimal  policy.  Our  main  intent 
with  Investigating  this  policy  rule  is  to  see  if  one  can  afford  to 
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ignore  the  sporadic  and  correlated  nature  of  the  warehouse  demand 
process.  An  evaluation  of  the  performance  of  this  policy  rule  is 
presented  in  Section  4. 

The  final  policy  rule  we  consider  is  the  Correlation-Adjusted 
Power  Approximation.  This  approximately  optimal  policy  rule,  which 
is  presented  in  Section  1,  takes  into  consideration  the  correlated 
nature  of  the  warehouse  demand  process.  An  evaluation  of  its  per- 
formance is  given  in  Section  5. 

3.2  Experimental  Design 

In  order  to  empirically  evaluate  the  performance  of  alternative 
(s,S)  policies  for  the  wholesale  warehouse  inventory  model,  we  need 
to  establish  an  experimental  design.  This  design  requires  the  speci- 
fication both  of  warehouse  parameters  and  of  store  parameters. 

3.2.1  Warehouse  Parameters 

We  consider  a ful 1 -factorial  72-item  warehouse  inventory  system. 

The  parameter  values  chosen  are  given  in  Table  3.2  . They  were  chosen 
to  be  consistent  with  those  used  in  previous  studies  by  Ehrhardt  (1976), 
Kaufman  (1977),  and  MacCormick(1974)  . The  four  values  for  mean  ware- 
house demand  are  4,  8,  12,  and  16  . The  variance-to-mean  ratio  of 
the  warehouse  demand  process  is  9 . Three  values,  L = 0,  2,  and  4 , 
are  assigned  to  leadtime.  Since  the  cost  function  is  linear  in  the 
parameters  K , h , and  p , the  value  of  unit  holding  cost  is  set 
at  unity.  Essential  changes  in  costs  arise  only  with  changes  in  the 
ratios  K/h  and  p/h  . The  penalty  cost  values  are  p = 4,  9,  and  99  . 
The  setup  cost  values  are  K = 32  and  64  . 


Warehouse  Parameters 


PARAMETER 

PARAMETER 

SETTINGS 

NUMBER 

OF 

SETTINGS 

Mean  Demand,  yw 

4,  8,  12,  16 

4 

Demand  Variance/  Mean, 

2. 

a / u 

w w 

9 

1 

Delivery  Leadtime,  L 

0,  2,  4 

3 

Unit  Holding  Cost,  h 

1 

1 

Unit  Penalty  Cost,  p 

4,  9,  99 

3 

Ordering  Setup  Cost,  K 

32,  64 

2 

3.2.2  Store  Parameters 

To  achieve  the  warehouse  demand  parameters  given  in  Table  3.1  , 
we  must  appropriately  specify  the  number  of  stores  and  the  demand  dis- 
tribution and  policy  parameters  at  each  of  the  stores.  In  addition, 
specification  of  the  store  parameters  will  induce  a warehouse  autocor- 
relation function.  To  aid  the  subsequent  analysis,  we  specify  the 
store  parameters  so  that  the  autocorrelation  function  of  the  warehouse 
demand  process  is  the  same  for  each  of  the  72-items  in  the  system. 


. . . i. 
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To  start,  we  consider  a simplified  environment,  specifically, 
one  in  which  each  store  employs  the  same  (s,S)  policy  and  has  the 
same  demand  distribution.  We  assume  that  each  store's  demand  distribu- 
tion is  a negative  binomial  distribution,  with  parameters  r and  p, 
where  r > 0 and  0 < p < 1 : 

(3.1)  4>(x)  = (l-p)rpXF(r+x)/[r(x+l  )r(r)]  for  x = 0,1,2,...  , 
yielding  a mean 

(3.2)  us  = rp/(l-p) 
a variance 

(3.3)  o2s  = rp/(l-p)2 
and  a variance-to-mean  ratio 

(3.4)  a2/ ps  = (1-p)"1 

It  is  easy  to  see  that  the  mean  and  variance  completely  characterize 
this  distribution. 

For  given  values  of  N and  D we  can  determine  values  for  y$ 

2 2 
and  a that  will  Induce  desired  values  of  y,  and  o . From  Lemma 
s w w 

2.1  and  Lemma  2.2  we  have 

(3.5)  Ps  = Ur  = Ww/N  * 

2 

The  problem  of  finding  the  appropriate  value  for  , however,  is 
much  more  difficult.  In  fact,  it  may  not  exist  at  all  if  D is  too 


. i. 
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large.  To  illustrate,  using  Lemma  2.1  and  Lemma  2.2  , we  can  write 
o o D 

(3.6)  o;  = (aVN)  - [2y  £ km(k)]/[N(l  + M(D))] 

S W w k=0 

? 

It  can  be  shown  that  o$  is  decreasing  in  D . For  the  negative 

2 

binomial  distribution  to  be  well-defined,  we  must  have  a$  > ys  . 

Upon  examination  of  (3.6)  we  see  that  it  is  entirely  possible  that 

a given  value  of  D can  be  too  large  for  this  inequality  to  hold. 

By  the  definitions  given  in  (2.3)  and  (2.4)  , the  renewal 

functions  m(*)  and  M(*)  depend  on  the  probability  mass  function 

$(•)  , or  equivalently,  for  the  negative  binomial  distribution  on 

Ps  and  . Thus  the  direct  solution  of  (3.6)  is  impossible. 

Nonetheless,  we  do  know  that  if  a solution  does  exist  it  must  be  con- 

tained  in  the  interval  (y  , o' /N)  , and  numerical  search  techniques 

s w 

can  be  employed  to  determine  an  approximate  solution  if  it  exists. 

Since  we  have  postulated  that  each  store  employs  the  same  (s,S) 
policy  and  has  the  same  demand  distribution,  we  have  from  Lemma  2.2 


(3.7) 


Pw(j)  = Pr(j) 


for  j = 1,2,3,...  . 


Thus  p (j)  is  independent  of  N for  all  lags  j . From  Lemma  2.1 
w 

2 

we  observe  that  p (j)  is,  therefore,  a function  of  us  , o$  , and 

2 

D . Thus  by  holding  y$  , o$  , and  D fixed,  we  can  change  the 

2 

values  of  y and  o proportionally  and,  at  the  same  time,  keep 
w w 

Pw ( j ) constant  for  all  j simply  by  changing  N in  the  appropriate 

2 

manner.  For  example,  to  double  the  values  of  yw  and  ow  and  keep 
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Pw(j)  fixed  for  all  j , we  simply  double  the  number  of  stores  N . 

In  this  investigation  we  consider  two  warehouse  demand  environ- 
ments. The  store  parameters  corresponding  to  these  demand  environ- 
ments are  given  in  Table  3.3  . Both  demand  environments  assume  a 
value  of  D = 8 and  a negative  binomial  demand  distribution  at  each 
store.  The  demand  environments  are  characterized  by  the  number  of 
stores  in  each,  with  one  having  four  times  as  many  stores  as  the 
other.  In  the  few-stores  environment  each  store  has  a mean  of  4 , 
a variance  of  12.80  , and  a replenishment  variance  of  36  , while 
for  the  many-stores  environment,  the  values  are,  respectively,  1 , 

1.70  , and  9 . Thus,  warehouse  mean  demand  values  of  4,  8,  12,  and 
16  correspond  to  N = 1,2,3,  and  4 for  the  few-stores  environment 
and  correspond  to  N = 4,8,12,  and  16  in  the  many-stores  environment. 


Table  3.3 
Store  Parameters 

(D  = 8 , negative  binomial  demand  distribution  at  each  store) 


DEMAND  ENVIRONMENT 

N 

ps 

2 

°s 

2 

°r 

Few  Stores 

V /4 
w 

4 

12.80 

36 

Many  Stores 

uw 

1 

1.70 

9 

The  warehouse  autocorrelation  functions  for  both  demand  environ- 
ments are  listed  in  Table  3.4  up  to  lag  four.  For  the  few-stores 
environment  the  autocorrelation  function  has  a value  of  -0.30  at 
lag  one  and  rapidly  approaches  zero  for  lags  greater  than  one.  For 
the  many-stores  environment  the  autocorrelation  function  has  a value 
of  -0.11  at  lag  one  and  slowly  goes  to  zero  as  the  lag  number 
increases.  The  difference  between  the  autocorrelation  functions  arises 
from  the  different  values  for  the  mean  and  variance  of  the  store's 
demand  distribution. 


Table  3.4 

Warehouse  Autocorrelation  Functions 


DEMAND  ENVIRONMENT 

p„0) 

PW(2) 

P„(3> 

P (4) 

Kwv  ' 

Few  Stores 

-0.30 

-0.05 

+0.03 

+0.01 

Many  Stores 

-0.11 

-0.10 

-0.09 

-0.07 

Figures  3.1  and  3.2  clearly  demonstrate  the  sporadic  nature 
of  the  warehouse  probability  mass  functions  for  the  few-stores  environ- 
ment. Note  for  N = 1 the  high  probability  of  zero  demand  and  the 
monotonical ly  decreasing  curve  for  demand  values  beyond  D+l  . When 
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N = 2 notice  a slight  increase  in  the  probability  mass  begins  at 
2(D+1)  . This  is  the  effect  of  the  second  store.  For  N = 3 and 
N = 4 we  still  observe  this  bulge  which  begins  at  2(D+1)  ; however, 
it  becomes  more  pronounced  when  N = 4 . The  effect  of  more  stores 
is  to  shift  more  mass  into  the  tail  of  the  distribution. 


3.2.3  The  Assumption  of  Identical  Stores 

In  Subsection  3.2.2  we  presented  an  experimental  design  in 
which  each  store  employs  the  same  (s,S)  policy  and  has  the  same 
demand  distribution.  One  may  suspect  that  the  warehouse  demand  pro- 
cess would  be  significantly  affected  by  this  assumption  of  identical 
stores.  We  have,  therefore,  computed  correlation  coefficients  and 
the  probability  of  zero  demand  for  the  warehouse  demand  process  under 
several  designs  with  non-identical  stores. 

Let  D.j  and  y..  denote,  respectively,  the  value  of  D and  the 
mean  of  the  demand  distribution  at  store  i . The  store  parameter 
values  corresponding  to  several  designs  with  non-identical  stores  are 
given  in  Table  3.5  . For  each  design  the  warehouse  demand  process 
has  a mean  of  8 and  a variance-to-mean  ratio  of  11  . In  addition, 

the  variance-to-mean  ratio  of  each  store's  replenishment  process  is 
also  set  at  11  . Note  that  these  variance-to-mean  ratios  have  been 
increased  from  the  previous  values  of  9 ; this  was  done  in  part  to 
allow  for  higher  D values.  For  each  demand  environment  there  are 
3 designs.  The  first  design  has  identical  stores,  the  second  has 
identical  D values  but  two  different  mean  demand  values,  and  the 
third  has  identical  mean  values  but  two  different  D values. 
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Table  3.5 

Several  Designs  With  Non-Identical  Stores 


DEMAND  ENVIRONMENT 

DESIGN 

VALUES  OF  D 

VALUES  OF  us 

Few  Stores 

I 

CO 

II 

C\J 

o 

II 

o 

M1  = U2  = 4 

II 

D1  = D2  = 8 

U1  = 2;  u2  = 6 

III 

D]  =6;  D2  = 10 

W]  = u2  = 4 

Many  Stores 

IV 

D1  = •••  = °8  = 8 

W1  = m8  = 1 

V 

o 

II 

II 

o 

00 

II 

00 

Pi  = •••  = P4  = 0.5; 

P5  = ...  = P8  = 1.5 

VI 

D-j  = ...  = D^  = 6; 

Pl  = ...  = U8  = ! 

D5  = ...  = Dg  = 10 

Table  3.6  gives  the  probability  that  warehouse  demand  in  period  i , 
Zi  , is  zero  along  with  several  autocorrelation  coefficients  for  each  of 
the  designs  given  in  Table  3.5  . 

Comparing  the  autocorrelation  coefficients  of  designs  I and  IV 
with  those  given  in  Table  3.4  , we  see  that  increasing  o /u  from  9 

W W 

to  11  results  in  less  correlation  between  warehouse  demands  but  similar 
behavior  of  the  autocorrelation  function.  By  comparing  design  I with 
designs  II  and  III  and  design  IV  with  designs  V and  VI  , we  see 
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that  no  significant  changes  have  occurred  either  in  the  probability 
of  zero  demand  or  in  the  autocorrelation  coefficients.  These  results 
suggest  that  the  assumption  of  identical  stores  is  not  a critical 
factor  in  influencing  the  nature  of  the  warehouse  demand  process. 


Table  3.6 

Probability  of  Zero  Demand  and  Autocorrelation 
Coefficients  for  the  Warehouse  Demand  Process 
Linder  Several  Non-Identical  Store  Designs 


DEMAND 

ENVIRONMENT 

DESIGN 

Pr[Z.  = 0] 

PW(D 

P„(2> 

P (3) 
pw'  ' 

P (4) 
w ' 

Few  Stores 

I 

0.488 

-0.06 

0.00 

+0.01 

II 

0.490 

-0.20 

-0.03 

0.00 

0.00 

III 

0.480 

-0.21 

-0.07 

+0.01 

+0.01 

Many  Stores 

n 

0.473 

-0.07 

-0.06 

-0.05 

-0.05 

0.475 

-0.09 

-0.07 

-0.05 

-0.03 

H 

0.465 

-0.07 

-0.06 

-0.06 

-0.05 

i 

3.3  Item  Operating  Characteristics 

In  Sections  4 and  5 the  evaluation  of  performance  of  alternative 
(s,S)  policy  rules,  under  the  experimental  design  described  above,  is 
made  by  comparing  their  operating  characteristic  values  with  those 
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obtained  using  best  (s,S)  policies  in  the  same  environment.  Operating 
characteristic  values  for  each  item  in  the  72-item  system  are  collected 
for  5000  periods.  The  operating  characteristics  we  consider  are  the 
expected  values  per  period  of  period-end  inventory,  backlog  quantity, 
frequency  of  period-end  backlogs,  replenishment  quentity,  frequency 
of  replenishment  and  the  total  cost  incurred. 


3.4  The  Multi -I tern  System 

Scientific  techniques  for  inventory  control  are  generally  applied 
to  systems  of  many  items.  Hence  in  Sections  4 and  5 we  also  com- 
bine the  results  from  two  72-item  systems  into  two  multi-item  systems. 
Since  management  generally  assesses  the  performance  of  control  techniques 
by  observing  indices  that  are  aggregate  operating  characteristics  [Wagner 
(1962)]  , certain  aggregate  characteristics  have  been  computed. 

The  operating  characteristics  of  the  multi-item  system  have  been 
measured  by  aggregating  the  sample  values  of  the  corresponding  charac- 
teristics for  each  item  in  the  system.  The  aggreqate  of  average  total 
cost  per  period  is  computed  as  the  arithmetic  sum  of  the  corresponding 
costs  for  each  item.  The  components  of  total  cost  for  inventory  storage, 
backlog  penalty,  and  ordering  replenishments  are  similarly  computed. 

The  aggregate  backlog  and  replenishment  frequencies  are  arithmetic  aver- 
ages of  the  corresponding  frequencies  observed  for  each  item  in  the 
system.  Since  the  unit  inventory  holding  cost  for  all  items  is  unity, 
the  average  number  of  units  in  inventory  at  period-end  is  numerically 
identical  to  the  aggregate  holding  cost  per  period.  Finally,  a weighted 
sum  of  the  average  quantity  backlogged  per  period  to  a weighted  sum  of 
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the  (exact)  mean  values  of  demand.  The  weights  used  in  both  the 
numerator  and  denominator  of  the  ratio  are  the  unit  cost  of  backlog- 


ging demand  for  the  respective  item. 


4.  EVALUATION  OF  A POLICY  RULE  THAT  IGNORES  SPORADIC  AND  CORRELATED 
DEMAND  PROPERTIES 

In  this  section  we  evaluate  the  performance  of  the  I.I.D.  optimal 
policy  introduced  in  Table  3.1  . Recall  that  under  this  policy  rule 
we  compute  optimal  policies  for  items  having  the  warehouse  cost  struc- 
ture, but  with  an  independent  and  identically  distributed  demand  pro- 
cess following  a negative  binomial  distribution  having  the  same  mean 
and  variance  as  the  warehouse  demand  process.  The  policies  are  com- 
puted by  the  algorithm  of  Veinott  and  Wagner  (1965)  as  programmed  by 
Kaufman  (1976)  . Essentially,  this  policy  rule  ignores  the  sporadic 
and  correlated  nature  of  the  warehouse  demand  environment.  The  pur- 
pose of  testing  the  I.I.D.  optimal  policy  rule  in  the  warehouse  demand 
environment  is  to  see  if  one  can  afford  to  ignore  these  properties. 

It  is  particularly  important  that  this  issue  be  studied,  since  practi- 
tioners often  make  this  kind  of  simplifying  assumption. 

We  find  that  the  use  of  the  I.I.D.  optimal  policy  rule  in  our 
warehouse  demand  environment  results  in  significant  degradation  in  per- 
formance. Specifically,  the  total  cost  performance  of  the  I.I.D.  optimal 
policies  is  about  8%  above  that  of  the  best  (s,S)  policies,  as 
defined  in  Table  3.1  . This  degradation  is  the  result  of  holding  too 
much  inventory. 

4.1  The  Performance  of  I.I.D.  Optimal  Policies 


Table  4.1  summarizes  the  performance  of  the  I.I.D.  optimal  policies 
in  both  warehouse  demand  environments.  The  components  of  aggregate 
average  total  cost  per  period  are  listed  and  are  compared  with  the 


Summary  of  the  Performance  of  the  I.I.D.  Optimal  Policies  for  Two 
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Numbers  in  parentheses  are  percentages  of  total  cost. 

Numbers  in  brackets  are  percentage  differences  in  cost  components  over  best  values 


corresponding  cost  when  the  systems  are  controlled  by  best  (s,S) 
policies.  Notice  that  the  total  cost  performance  of  the  I.I.D.  optimal 
policies  is  8.8%  and  7.1%  higher  than  the  best  values  for  the  few- 
stores  environment  and  the  many-stores  environment,  respectively.  The 
cause  of  this  degradation  in  performance  is  clearly  the  result  of  hold- 
ing too  much  inventory.  Specifically,  31%  more  inventory  is  held  in 
the  few-stores  environment,  and  25%  more  in  the  many-stores  environ- 
ment. In  addition,  inventory  costs  account  for  approximately  two- 
thirds  of  the  total  cost  in  both  environments.  Observe  that  backlog 
and  replenishment  costs  are  significantly  lower  in  both  environments; 
however,  the  savings  over  best  values  are  not  subtantial  enough  to 
offset  the  rise  in  inventory  holding  costs. 

Insight  into  the  cost  component  differences  is  gained  by  compar- 
ing the  policy  parameter  values  s and  S of  the  best  policies  with 
those  of  the  I.I.D.  optimal  policies.  Table  4.2  lists  policy  param- 
eter values  for  a subset  of  the  72-item  few-stores  environment.  Note 
that  for  the  best  policies  D = S-s  is  smaller  for  all  items  except 
one,  s is  smaller  for  all  items  except  those  with  zero  leadtime 
and  small  penalty  cost,  and  S is  smaller  with  the  exception  of  two 
items.  This  subset  is  representative  of  all  72-items  and  both  environ- 
ments. Thus,  in  general,  D and  S tend  to  be  smaller,  and  s is 
smaller  except  for  those  items  with  zero  leadtime  and  small  unit  penalty 
costs.  A complete  comparison  of  the  policy  parameters  is  given  in 
Technical  Report  15.  The  smaller  D values  account  for  the  larger 
replenishment  costs  for  the  best  policies,  and  the  smaller  s values 
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for  those  items  with  high  penalty  costs  result  in  increased  aggregate 
backlog  costs.  However,  the  S values  which  are  also  smaller,  par- 
ticularly so  for  those  items  with  a larger  penalty  cost  and  positive 
leadtime,  result  in  considerable  holding  cost  savings. 

Detailed  comparisons  of  the  systems  under  the  two  decisions  rules, 
expanding  on  Table  4.1  , are  found  in  Tables  4.3  and  4.4  showing 
percentages  above  best  values  for  each  cost  component  and  parameter 
setting  using  the  I.I.D.  optimal  policies.  The  total  cost  data 
reveal  that  the  I.I.D.  optimal  policies  are  close  to  best  (s,S) 
policies  only  for  those  items  with  small  penalty  costs  or  those  with 
zero  leadtime.  High  penalty  costs  and  large  leadtimes  are  particularly 
troublesome,  and  to  a somewhat  lesser  degree,  so  are  low  means  and 
small  setup  costs.  Inventory  costs  are  considerably  above  best  values 
for  all  parameter  settings,  especially  for  a high  penalty  cost  or  a 
large  leadtime.  Backlog  costs  are  below  best  values  for  nearly  all 
of  the  parameter  settings.  Note,  in  particular,  that  backlog  costs 
are  considerably  less  when  the  penalty  cost  is  99  or  when  the 
leadtime  is  4 . Replenishment  costs  are  below  best  values  for  all 
parameter  settings  and,  similar  to  backlog  costs,  have  the  largest 
discrepancy  for  those  items  with  high  penalty  costs  or  large  leadtimes. 
Cost  component  differences  are  similar  for  both  environments,  the 
many-stores  environment  exhibiting  slightly  smaller  differences. 

The  percentage  apportionment  of  aggregate  costs  per  period  for 
various  parameter  classifications  is  shown  in  Tables  4.5  and  4.6 
for  best  (s,S)  policies  and  for  the  I.I.D.  optimal  (s,S)  policies. 
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Table  4.4 
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Inventory  costs  for  the  I.I.D.  optimal  policies  increase  more  dramati- 
cally for  higher  penalty  costs  and  for  positive  leadtimes  than  for 
best  policies.  Backlog  costs,  on  the  other  hand,  decrease  far  more 
sharply  as  penalty  costs  increase.  Notice  that  as  the  leadtime 
increases,  the  backlog  costs  increase  for  the  best  policies  and 
decrease  for  the  I.I.D.  optimal  policies. 

Tables  4.7  and  4.8  show  the  values  of  other  operating  charac- 
teristics of  the  systems  under  best  policies  and  I.I.D.  optimal  poli- 
cies. Backlog  and  replenishment  frequency  are  higher  for  the  best 
policies.  Backlog  frequency  increases  only  slightly  in  leadtime  for 
best  policies  and  decreases  sharply  for  I.I.D.  optimal  policies. 

Higher  setup  costs  lead  to  less  frequent  backlogs  for  the  best  poli- 
cies and  more  frequent  backlogs  for  the  I.I.D.  optimal  policies. 
Percentages  above  best  values  for  the  operating  characteristics  of 
the  I.I.D.  optimal  policies  are  given  in  Tables  4.9  and  4.10  . The 
results  are  consistent  with  trends  noted  for  the  cost  components. 

The  above  performance  summary  indicates  that  the  I.I.D.  optimal 
policy  rule  is  unsatisfactory,  especially  for  those  items  with  a high 
penalty  cost  or  a long  leadtime.  In  general , the  I.I.D.  optimal  poli- 
cies hold  too  much  inventory  and  order  too  infrequently.  Hence,  it 
can  be  very  costly  to  ignore  the  sporadic  and  correlated  nature  of  the 
warehouse  demand  process.  What  is  needed  is  a policy  rule  that  can 
adjust  itself  to  these  demand  characteristics  in  a simple  way.  In 
Section  5 we  present  a policy  rule  which  adjusts  for  correlation  in 
a demand  process. 
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5.  THE  POWER  APPROXIMATION  IN  A WAREHOUSE  DEMAND  ENVIRONMENT 


In  the  previous  section  we  observed  that  the  correlated  and 
sporadic  nature  of  the  warehouse  demand  process  cannot  be  overlooked 
without  serious  degradation  in  system  performance.  Klincewicz  (1976) 
studied  the  effectiveness  of  the  Power  Approximation  in  a particular 
sporadic  demand  environment  in  which  demands  are  independent  and 
identically  distributed.  He  studied  demand  distributions  that  have 
considerable  probability  of  no  demand  in  a period  and  also  have  a 
high  variance-to-mean  ratio.  His  results  indicate  that  the  Power 
Approximation  performance  in  such  an  environment  is  close  to  optimal. 
Although  more  pronounced,  the  sporadic  effect  in  the  warehouse  demand 
environment  is  similar.  For  this  reason  we  have  concentrated  on  the 
correlation  effects  of  the  warehouse  demand  process  on  the  policy 
parameters. 

In  Subsection  5.1  we  present  an  adaptation  of  the  Power  Approxi- 
mation to  a correlated  demand  environment.  The  modified  algorithm 
requires  only  knowledge  of  the  mean,  variance,  and  autocorrelation 
function  up  to  lag  L . In  Subsection  5.2  we  evaluate  the  performance 
of  the  Correlation-Adjusted  Power  Approximation.  We  find  its  perform- 
ance to  be  close  to  the  performance  of  the  best  (s,S)  policies. 
Specifically,  for  the  few-stores  and  many-stores  environments  the  poli- 
cies yield  costs  of  only  2.6%  and  3.4%  , respectively,  above 
expected  costs  for  the  best  policies. 


56 


r 


5 . 1 Adaptation  of  the  Power  Approximation  to  Correlated  Demands 

The  Power  Approximation  is  an  algorithm  for  computing  approxi- 
mately optimal  values  for  (s,S)  using  only  the  mean  and  variance 
of  demand.  Ehrhardt  (1976)  has  demonstrated  that  the  algorithm  is 
a highly  accurate  approximation  to  optimal  (s,S)  policies.  For 
this  reason  the  performance  of  Power  Approximation  policies  in  a 
warehouse  demand  environment  will  closely  resemble  that  of  the  I.I.D. 
optimal  policies.  Nonetheless,  we  can  modify  the  Power  Approximation 
for  use  in  a correlated  demand  environment. 

In  the  Power  Approximation  expressions  for  s and  D we  find 
the  term 


(5.1) 


= (L  + 1 )o2 


2 

where  a represents  the  variance  of  the  dem&nd'process . Since  the 
Power  Approximation  was  derived  for*"7ndependent,  identically  distributed 
demands,  the  term  represents  the  variance  of  demand  over  L+l 

periods.  We  modify  the  Power  Approximation  by  simply  replacing  (5.1) 
with  the  more  general  expression 


(5.2) 


JL+1 


= [(L+l) 


+ 2 l 
j=l 


(L+l-j)p(j)]o 


where  p(j)  is  the  autocorrelation  of  the  process  at  lag  j . We 
refer  to  the  Power  Approximation  with  (5.2)  replacing  (5.1)  as 


i 


the  Correlation-Adjusted  Power  Approximation. 

Upon  examination  of  the  autocorrelation  functions  for  both  demand 
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environemnts  in  Table  3.2  , it  is  clear  that  the  correlation-adjusted 

2 

variance  over  L+l  periods  will  be  less  than  (L+l)o  . Table  5.1 

2 

lists  this  adjusted  variance  for  several  values  of  a and  L and 
for  both  demand  environments.  Notice  the  significant  reduction  in 
the  variance  over  L+l  periods,  especially  for  the  few-stores  environ- 
ment, where  reductions  are  in  the  neighborhood  of  40  to  50  percent. 
Table  5.2  compares  Power  Approximation  policies  using  (5.1) 

and  (5.2)  for  a subset  of  the  72-item  systems.  The  values  of  D 

2 

are  only  1 or  2 units  smaller,  since  is  raised  to  a very 

small  power  in  the  expression  for  D . However,  s has  decreased 
substantially,  especially  when  the  penalty  cost  is  99  . By  comparing 
with  best  parameter  values,  we  see  that  for  most  items  the  changes  in 
D and  s are  in  the  appropriate  direction.  Note  also  that  the  poli- 
cies change  more  dramatically  for  the  few-stores  environment. 

5.2  The  Performance  of  Correlation-Adjusted  Power  Approximation  Policies 

Table  5.3  summarizes  the  performance  of  the  Correlation-Adjusted 
Power  Approximation  policies  for  both  warehouse  demand  environments. 

The  components  of  aggregate  average  total  cost  per  period  are  listed 
and  are  compared  with  the  corresponding  costs  when  the  systems  are  con- 
trolled by  best  (s,S)  policies.  Notice  the  total  cost  component  is 
only  2.6%  above  best  (s,S)  policy  performance  for  the  few-stores 
environment,  and  3.4%  above  for  the  many-stores  environment.  Inventory 
holding  costs  are  above  best  values,  while  backlog  and  replenishment 
costs  are  below  best  values.  Recall  that  the  I.I.D.  optimal  policies 
displayed  the  same  characteristics;  however,  here  the  differences  are 


Uncorrelated  and  Correlated  Demands 


258 


Best 


135 


61 


less  severe.  Substantial  improvement  has  been  made  in  the  inventory 
cost  component  due  to  lower  values  of  S . In  particular,  inventory 
costs  for  the  few-stores  environment  are  down  from  31%  to  10.5% 
above  best  values.  Replenishment  costs  for  both  environments  are, 
respectively,  11.5%  and  12.5%  below  best  values.  Since  D is 
virtually  unchanged,  these  percentages  are  only  slightly  smaller  than 
the  corresponding  percentages  when  using  the  I.I.D.  optimal  policies. 
Backlog  costs  are,  respectively,  1%  and  5.5%  below  best  values 
for  the  few-stores  environment  and  the  many-stores  environment,  and 
are  much  closer  to  best  values  than  the  corresponding  backlog  costs 
using  the  I.I.D.  optimal  policies.  This  can  be  attributed  to  the 
lower  values  of  s . 

Tables  5.4  and  5.5  show  percentages  above  best  values  for 
each  cost  component  and  parameter  setting  suing  the  Correlation- 
Adjusted  Power  Approximation  policies.  The  total  cost  component 
reveals  that,  with  the  exception  of  those  items  with  high  penalty 
cost,  the  Correlation-Adjusted  Power  Approximation  performance  is 
close  to  best  performance.  Comparing  with  Tables  4.3  and  4.4  , 
we  note  that  the  severe  degradation  for  large  leadtimes  with  the  I.I.D. 
optimal  policies  is  not  longer  present.  That  is,  our  revised  calcula- 
tion of  leadtime  variance  has  effected  a substantial  improvement  for 
those  items  with  positive  leadtime. 

The  percentage  apportionment  of  aggregate  costs  per  period  for 
various  parameter  classifications  is  shown  in  Tables  5.6  and  5.7 
for  best  (s,S)  policies  and  for  Correlation-Adjusted  Power  Approxi- 
mation policies.  Note  that  backlog  costs  increase  slightly  as 


Percentage  Above  Best  Costs  Per  Period  for  a 72-Item  Warehouse  Inventory  System 


32.3  36.5  16.5 
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leadtime  increases,  which  is  now  in  agreement  with  best  policy 
performance. 

Tables  5.8  and  5.9  show  the  values  of  other  operating 
characteristics  of  the  systems  under  best  policies  and  Correlation- 
Adjusted  Power  Approximation  policies.  Note  that  for  the  Correlati 
Adjusted  Power  Approximation,  backlog  frequency  is  slightly  depen- 
dent on  setup  costs,  whereas  best  policies  exhibit  no  such  depen- 
dence. Percentages  above  best  values  for  the  operating  characters 
tics  of  the  Correlation-Adjusted  Power  Approximation  policies  are 
given  in  Tables  5.10  and  5.11  . Note  backlog  frequencies  are 
above  best  values  except  for  high  penalty  cost  items,  in  which 
backlog  frequency  is  substantially  below  the  best  value. 
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6.  CONCLUSIONS 


In  this  section  we  summarize  the  principal  results  of  this 
report  and  outline  topics  to  be  investigated  further. 

In  Section  2 we  presented  a detailed  analysis  of  the  demand 
process  for  our  wholesale  warehouse  inventory  model.  We  were  able 
to  derive  useful  expressions  for  its  mean,  variance,  and  autocorre- 
lation function  in  terms  of  the  policy  parameters  and  demand  distri- 
bution at  each  store. 

We  considered  several  (s,S)  policy  rules  for  the  wholesale 
warehouse.  In  Section  4 we  examined  a policy  rule  which  ignores 
both  the  sporadic  and  the  correlated  nature  of  the  warehouse  demand 
process.  This  policy  rule  is  optimal  if  demands  are  independent 
and  identically  distributed  negative  binomial  random  variables  with 
a mean  and  variance  equal  to  that  of  the  warehouse  demand  process. 

We  found  the  performance  of  this  policy  rule  to  be  unsatisfactory 
in  several  warehouse  demand  environments;  particularly,  when  the 
penalty  cost  is  high  or  the  leadtime  is  large.  In  general,  the  policy 
rule  holds  too  much  inventory  and  orders  too  infrequently.  In  Sec- 
tion 5 we  examined  a policy  rule  which  takes  into  account  the  corre- 
lated nature  of  the  demand  process.  This  policy  rule  is  a simple 
modification  of  the  Ehrhardt  Power  Approximation.  The  only  demand 

information  required  is  the  mean  and  variance  of  demand  and  the 

2 

variance  of  demand  over  one  leadtime,  o^+-|  . We  found  that  this 
policy  performed  quite  well  in  the  same  warehouse  demand  environments. 

When  using  the  Correlation-Adjusted  Power  Approximation  in  this 
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study,  we  assumed  that  y,  a , and  o^+-|  are  known.  How  this 

policy  rule  performs  when  these  demand  parameters  must  be  estimated 

from  a limited  demand  history  is  an  important  issue  still  to  be 

investigated.  The  main  difficulty  foreseen  is  the  accurate  estima- 
2 

tion  of  o^+-|  . One  possible  method  is  to  estimate  the  autocorrela- 
tion function  up  to  lag  L , then  let 

(6.1)  , = [(L+l)  + 2 l (L+l-j )p(j )]o^  . 

j = l 

If  the  demand  history  is  limited,  however,  then  the  estimate  p(j) 
is  likely  to  be  very  unstable  for  higher  lags.  One  possible  solu- 
tion is  to  consider  estimating  only  the  first-order  autocorrelation 
and  setting  the  others  equal  to  zero.  As  a preliminary  step,  we 
have  investigated  the  performance  of  the  Correlation-Adjusted  Power 
Approximation  when  we  know  the  first-order  autocorrelation  p(l)  , 
and  substitute  p(j)  = 0 for  all  j > 1 in  expression  (5.2)  . A 
summary  of  its  performance  is  provided  in  Table  6.1  . If  we  compare 
these  results  with  those  given  in  Table  5.3  , using  all  of  the 
higher-order  autocorrelations , we  see  that  only  a slight  degradation 
in  performance  has  occurred.  hese  results  are  encouraging  and  suggest 
we  may  indeed  only  need  to  estimate  p(l)  . 

Finally,  the  robustness  of  the  Correlation-Adjusted  Power  Approxi- 
mation in  other  types  of  correlated  demand  environments  needs  to  be 
studied.  For  example,  its  performance  in  an  environment  where  demand 
follows  a Markov  process  will  be  investigated. 
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NOTE:  Numbers  in  parentheses  are  percentages  of  total  cost 
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